Identification and Optimization of RNA-Splicing Modulators as Huntingtin Protein-Lowering Agents for the Treatment of Huntington's Disease.

Huntington's disease (HD) is caused by an expanded CAG trinucleotide repeat in exon 1 of the huntingtin (HTT) gene. We report the design of a series of HTT pre-mRNA splicing modulators that lower huntingtin (HTT) protein, including the toxic mutant huntingtin (mHTT), by promoting insertion of a pseudoexon containing a premature termination codon at the exon 49-50 junction. The resulting transcript undergoes nonsense-mediated decay, leading to a reduction of HTT mRNA transcripts and protein levels. The starting benzamide core was modified to pyrazine amide and further optimized to give a potent, CNS-penetrant, and orally bioavailable HTT-splicing modulator 27. This compound reduced canonical splicing of the HTT RNA exon 49-50 and demonstrated significant HTT-lowering in both human HD stem cells and mouse BACHD models. Compound 27 is a structurally diverse HTT-splicing modulator that may help understand the mechanism of adverse effects such as peripheral neuropathy associated with branaplam.

Lactase deficiency in Russia: multiethnic genetic study.

BACKGROUND: Lactase persistence-the ability to digest lactose through adulthood-is closely related to evolutionary adaptations and has affected many populations since the beginning of cattle breeding. Nevertheless, the contrast initial phenotype, lactase non-persistence or adult lactase deficiency, is still observed in large numbers of people worldwide. METHODS: We performed a multiethnic genetic study of lactase deficiency on 24,439 people, the largest in Russia to date. The percent of each population group was estimated according to the local ancestry inference results. Additionally, we calculated frequencies of rs4988235 GG genotype in Russian regions using the information of current location and birthplace data from the client's questionnaire. RESULTS: The attained results show that among all studied population groups, the frequency of GG genotype in rs4988235 is higher than the average in the European populations. In particular, the prevalence of lactase deficiency genotype in the East Slavs group was 42.8% (95% CI: 42.1-43.4%). We also investigated the regional prevalence of lactase deficiency based on the current place of residence. CONCLUSIONS: Our study emphasizes the significance of genetic testing for diagnostics, i.e., specifically for lactose intolerance parameter, as well as the scale of the problem of lactase deficiency in Russia which needs to be addressed by the healthcare and food sectors.

Association of common genetic variants with body mass index in Russian population.

BACKGROUND: Overweight is the scourge of modern society and a major risk factor for many diseases. For this reason, understanding the genetic component predisposing to high body mass index (BMI) seems to be an important task along with preventive measures aimed at improving eating behavior and increasing physical activity. METHODS: We analyzed genetic data of a European cohort (n = 21,080, 47.25% women, East Slavs ancestry >80%) for 5 frequently found genes in the context of association with obesity: IPX3 (rs3751723), MC4R (rs17782313), TMEM18 (rs6548238), PPARG (rs1801282) and FTO (rs9939609). RESULTS: Our study revealed significant associations of FTO (rs9939609) (ß = 0.37 (kg/m2)/allele, p = <2 × 10-16), MC4R (rs17782313) (ß = 0.28 (kg/m2)/allele, p = 5.79 × 10-9), TMEM18 (rs6548238) (ß = 0.29 (kg/m2)/allele, p = 2.43 × 10-8) with BMI and risk of obesity. CONCLUSIONS: The results confirm the contribution of FTO, M4CR, and TMEM18 genes to the mechanism of body weight regulation and control.

Novel mutation in the MPZ gene causes early-onset but slow-progressive Charcot-Marie-Tooth disease in a Russian family: a case report.

Charcot-Marie-Tooth disease (CMT) is a genetically heterogeneous group of peripheral neuropathies most of which are associated with mutations in four genes including peripheral myelin protein-22 (PMP22), myelin protein zero (MPZ), gap junction protein beta1 (GJB1) and mitofusin2 (MFN2). This current case report describes the clinical and genetic characteristics of a 6-year-old male proband. A physical examination revealed muscular hypotonia. He started walking on his own at 18 months. A nerve conduction study with needle electromyography revealed conduction block. A novel MPZ mutation (c.398C > T, p.Pro133Leu) was revealed in the proband. This mutation was also found in the 32-year-old father of the proband. The father had had deformity of the feet and distal muscle weakness since childhood. The novel p.Pro133Leu pathogenic mutation was responsible for early onset but slowly progressive CMT1B. We assume that this site is an intolerant to change region in the MPZ gene. This variant in the MPZ gene is an important contributor to hereditary neuropathy with reduced nerve conduction velocity in the Russian population. This case highlights the importance of whole exome sequencing for a proper clinical diagnosis of CMT associated with a mutation in the MPZ gene.

GWAS of depression in 4,520 individuals from the Russian population highlights the role of MAGI2 (S-SCAM) in the gut-brain axis.

We present the results of the depression Genome-wide association studies study performed on a cohort of Russian-descent individuals, which identified a novel association at chromosome 7q21 locus. Gene prioritization analysis based on already known depression risk genes indicated MAGI2 (S-SCAM) as the most probable gene from the locus and potential susceptibility gene for the disease. Brain and gut expression patterns were the main features highlighting functional relatedness of MAGI2 to the previously known depression risk genes. Local genetic covariance analysis, analysis of gene expression, provided initial suggestive evidence of hospital anxiety and depression scale and diagnostic and statistical manual of mental disorders scales having a different relationship with gut-brain axis disturbance. It should be noted, that while several independent methods successfully in silico validate the role of MAGI2, we were unable to replicate genetic association for the leading variant in the MAGI2 locus, therefore the role of rs521851 in depression should be interpreted with caution.

Induction of Premature Cell Senescence Stimulated by High Doses of Antioxidants Is Mediated by Endoplasmic Reticulum Stress.

In our previous study, we found that high doses of several substances with antioxidant capacities (Tempol, resveratrol, diphenyleneiodonium) can cause genotoxic stress and induce premature senescence in the human mesenchymal stem cells (MSCs). Here, using whole-transcriptome analysis, we revealed the signs of endoplasmic reticulum stress and unfolded protein response (UPR) in MSCs stressed with Tempol and resveratrol. In addition, we found the upregulation of genes, coding the UPR downstream target APC/C, and E3 ubiquitin ligase that regulate the stability of cell cycle proteins. We performed the molecular analysis, which further confirmed the untimely degradation of APC/C targets (cyclin A, geminin, and Emi1) in MSCs treated with antioxidants. Human fibroblasts responded to antioxidant applications similarly. We conclude that endoplasmic reticulum stress and impaired DNA synthesis regulation can be considered as potential triggers of cell damage and premature senescence stimulated by high-dose antioxidant treatments.

Two novel PCDH19 mutations in Russian patients with epilepsy with intellectual disability limited to females: a case report.

BACKGROUND: Epilepsy with intellectual disability limited to females (Epileptic encephalopathy, early infantile, 9; EIEE9) is a rare early infantile epileptic encephalopathy characterized by an unusual X-linked inheritance: females with heterozygous mutations are affected, while hemizygous males are not. CASE PRESENTATION: We describe the clinical and molecular characteristics of 2 Russian patients with EIEE9 (females, ages 3 years and 7 years). In these patients seizures developed at the age of 3 years. Additionally, for our patients and for cases described in the literature we searched for a possible relationship between the type and localization of the mutation and the EIEE9 clinical phenotype. CONCLUSIONS: We identified two novel PCDH19 mutations in EIEE9 patients: a missense mutation in exon 1 (c.1236C > A, p.Asp412Glu) and a frameshift in exon 3 (c.2386_2387insGTCT, p.Thr796fs). We conclude that the age of seizure onset and the presence of intellectual disability may depend not on the type and localization of PCDH19 mutations, but on the X-inactivation status. The study also highlights the need to screen for EIEE9 among young female epilepsy patients.

Correction to Identification of Quinoline-Based RIP2 Kinase Inhibitors with an Improved Therapeutic Index to the hERG Ion Channel.

[This corrects the article DOI: 10.1021/acsmedchemlett.8b00344.].

Neuronal ceroid lipofuscinosis in the Russian population: Two novel mutations and the prevalence of heterozygous carriers.

BACKGROUND: Neuronal ceroid lipofuscinoses (NCLs) are a group of neurodegenerative disorders characterized by an accumulation of lipofuscin in the body's tissues. NCLs are associated with variable age of onset and progressive symptoms including seizures, psychomotor decline, and loss of vision. METHODS: We describe the clinical and molecular characteristics of four Russian patients with NCL (one female and three males, with ages ranging from 4 to 5 years). The clinical features of these patients include cognitive and motor deterioration, seizures, stereotypies, and magnetic resonance imaging signs of brain atrophy. Exome sequencing was performed to identify the genetic variants of patients with NCL. Additionally, we tested 6,396 healthy Russians for NCL alleles. RESULTS: We identified five distinct mutations in four NCL-associated genes of which two mutations are novel. These include a novel homozygous frameshift mutation in the CLN6 gene, a compound heterozygous missense mutation in the KCTD7 gene, and previously known mutations in KCTD7, TPP1, and MFSD8 genes. Furthermore, we estimated the Russian population carrier frequency of pathogenic and likely pathogenic variants in 13 genes associated with different types of NCL. CONCLUSION: Our study expands the spectrum of mutations in lipofuscinosis. This is the first study to describe the molecular basis of NCLs in Russia and has profound and numerous clinical implications for diagnosis, genetic counseling, genotype-phenotype correlations, and prognosis.

Identification of Quinoline-Based RIP2 Kinase Inhibitors with an Improved Therapeutic Index to the hERG Ion Channel.

RIP2 kinase was recently identified as a therapeutic target for a variety of autoimmune diseases. We have reported previously a selective 4-aminoquinoline-based RIP2 inhibitor GSK583 and demonstrated its effectiveness in blocking downstream NOD2 signaling in cellular models, rodent in vivo models, and human ex vivo disease models. While this tool compound was valuable in validating the biological pathway, it suffered from activity at the hERG ion channel and a poor PK/PD profile thereby limiting progression of this analog. Herein, we detail our efforts to improve both this off-target liability as well as the PK/PD profile of this series of inhibitors through modulation of lipophilicity and strengthening hinge binding ability. These efforts have led to inhibitor 7, which possesses high binding affinity for the ATP pocket of RIP2 (IC50 = 1 nM) and inhibition of downstream cytokine production in human whole blood (IC50 = 10 nM) with reduced hERG activity (14 µM).

Biopharmaceutical Informatics: supporting biologic drug development via molecular modelling and informatics.

OBJECTIVES: The purpose of this article is to introduce an emerging field called 'Biopharmaceutical Informatics'. It describes how tools from Information technology and Molecular Biophysics can be adapted, developed and gainfully employed in discovery and development of biologic drugs. KEY FINDINGS: The findings described here are based on literature surveys and the authors' collective experiences in the field of biologic drug product development. A strategic framework to forecast early the hurdles faced during drug product development is weaved together and elucidated using chemical degradation as an example. Efficiency of translating biologic drug discoveries into drug products can be significantly improved by combining learnings from experimental biophysical and analytical data on the drug candidates with molecular properties computed from their sequences and structures via molecular modeling and simulations. SUMMARY: Biopharmaceutical Informatics seeks to promote applications of computational tools towards discovery and development of biologic drugs. When fully implemented, industry-wide, it will enable rapid materials-free developability assessments of biologic drug candidates at early stages as well as streamline drug product development activities such as commercial scale production, purification, formulation, analytical characterization, safety and in vivo performance.

Quantifying the Risks of Asparagine Deamidation and Aspartate Isomerization in Biopharmaceuticals by Computing Reaction Free-Energy Surfaces.

Early identification of asparagine deamidation and aspartate isomerization degradation sites can facilitate the successful development of biopharmaceuticals. Several knowledge-based models have been proposed to assess these degradation risks. In this study, we propose a physics-based approach to identify the degradation sites on the basis of the free-energy barriers along the prechemical conformational step and the chemical reaction pathway. These contributions are estimated from classical and quantum mechanics/molecular mechanics molecular dynamics simulations. The computed barriers are compared to those for reference reactions in water within GNG and GDG sequence motifs in peptides (which demonstrate the highest degradation rates). Two major factors decreasing the degradation rates relative to the reference reactions are steric hindrance toward accessing reactive conformations and replacement of water by less polar side chains in the solvation shell of transition states. Among the potential degradation sites in the complementarity-determining region of trastuzumab and between two DK sites in glial cell-derived neurotropic factor, this method identified N30T, N55G, D102G, and D95K, respectively, in agreement with experiments. This approach can be incorporated in early computational screening of chemical degradation sites in biopharmaceuticals.

Computing the Free Energy Barriers for Less by Sampling with a Coarse Reference Potential while Retaining Accuracy of the Target Fine Model.

Proposed in this contribution is a protocol for calculating fine-physics (e.g., ab initio QM/MM) free-energy surfaces at a high level of accuracy locally (e.g., only at reactants and at the transition state for computing the activation barrier) from targeted fine-physics sampling and extensive exploratory coarse-physics sampling. The full free-energy surface is still computed but at a lower level of accuracy from coarse-physics sampling. The method is analytically derived in terms of the umbrella sampling and the free-energy perturbation methods which are combined with the thermodynamic cycle and the targeted sampling strategy of the paradynamics approach. The algorithm starts by computing low-accuracy fine-physics free-energy surfaces from the coarse-physics sampling in order to identify the reaction path and to select regions for targeted sampling. Thus, the algorithm does not rely on the coarse-physics minimum free-energy reaction path. Next, segments of high-accuracy free-energy surface are computed locally at selected regions from the targeted fine-physics sampling and are positioned relative to the coarse-physics free-energy shifts. The positioning is done by averaging the free-energy perturbations computed with multistep linear response approximation method. This method is analytically shown to provide results of the thermodynamic integration and the free-energy interpolation methods, while being extremely simple in implementation. Incorporating the metadynamics sampling to the algorithm is also briefly outlined. The application is demonstrated by calculating the B3LYP//6-31G*/MM free-energy barrier for an enzymatic reaction using a semiempirical PM6/MM reference potential. These modifications allow computing the activation free energies at a significantly reduced computational cost but at the same level of accuracy compared to computing full potential of mean force.

Quantitative exploration of the molecular origin of the activation of GTPase.

GTPases play a major role in cellular processes, and gaining quantitative understanding of their activation demands reliable free energy surfaces of the relevant mechanistic paths in solution, as well as the interpolation of this information to GTPases. Recently, we generated ab initio quantum mechanical/molecular mechanical free energy surfaces for the hydrolysis of phosphate monoesters in solution, establishing quantitatively that the barrier for the reactions with a proton transfer (PT) step from a single attacking water (1 W) is higher than the one where the PT is assisted by a second water (2 W). The implication of this finding on the activation of GTPases is quantified here, by using the ab initio solution surfaces to calibrate empirical valence bond surfaces and then exploring the origin of the activation effect. It is found that, although the 2 W PT path is a new element, this step is not rate determining, and the catalytic effect is actually due to the electrostatic stabilization of the pre-PT transition state and the subsequent plateau. Thus, the electrostatic catalytic effect found in our previous studies of the Ras GTPase activating protein (RasGAP) and the elongation factor-Tu (EF-Tu) with a 1 W mechanism is still valid for the 2 W path. Furthermore, as found before, the corresponding activation appears to involve a major allosteric effect. Overall, we believe that our finding is general to both GTPases and ATPases. In addition to the biologically relevant finding, we also provide a critical discussion of the requirements from reliable surfaces for enzymatic reactions.

Quantifying the mechanism of phosphate monoester hydrolysis in aqueous solution by evaluating the relevant ab initio QM/MM free-energy surfaces.

Understanding the nature of the free-energy surfaces for phosphate hydrolysis is a prerequisite for understanding the corresponding key chemical reactions in biology. Here, the challenge has been to move to careful ab initio QM/MM (QM(ai)/MM) free-energy calculations, where obtaining converging results is very demanding and computationally expensive. This work describes such calculations, focusing on the free-energy surface for the hydrolysis of phosphate monoesters, paying special attention to the comparison between the one water (1W) and two water (2W) paths for the proton-transfer (PT) step. This issue has been explored before by energy minimization with implicit solvent models and by nonsystematic QM/MM energy minimization, as well as by nonsystematic free-energy mapping. However, no study has provided the needed reliable 2D (3D) surfaces that are necessary for reaching concrete conclusions. Here we report a systematic evaluation of the 2D (3D) free-energy maps for several relevant systems, comparing the results of QM(ai)/MM and QM(ai)/implicit solvent surfaces, and provide an advanced description of the relevant energetics. It is found that the 1W path for the hydrolysis of the methyl diphosphate (MDP) trianion is 6-9 kcal/mol higher than that the 2W path. This difference becomes slightly larger in the presence of the Mg(2+) ion because this ion reduces the pKa of the conjugated acid form of the phosphate oxygen that accepts the proton. Interestingly, the BLYP approach (which has been used extensively in some studies) gives a much smaller difference between the 1W and 2W activation barriers. At any rate, it is worth pointing out that the 2W transition state for the PT is not much higher that the common plateau that serves as the starting point of both the 1W and 2W PT paths. Thus, the calculated catalytic effects of proteins based on the 2W PT mechanistic model are not expected to be different from the catalytic effects predicted using the 1W PT mechanistic model, which was calibrated on the observed barrier in solution and in which the TS charge distribution was similar to the that of the plateau (as was done in all of our previous EVB studies).

Addressing open questions about phosphate hydrolysis pathways by careful free energy mapping.

The nature and mechanism of phosphate hydrolysis reactions are of great interest in view of the crucial role of these reactions in key biological processes. Although it is becoming clearer that the ultimate way of resolving mechanistic controversies must involve reliable theoretical studies, it is not widely realized that such studies cannot be performed at present by using most existing automated ways and that only careful systematic studies can lead to meaningful conclusions. The present work clarifies the above point by considering the hydrolysis of phosphate monoesters. The clarification starts by defining the actual issues that should be addressed in careful studies and by highlighting the problems with studies that ignore the need for unique mechanistic definitions (e.g., works that confuse associative and dissociative pathways). We then focus on the analysis of the proton transfer (PT) pathways in phosphate hydrolysis and on recent suggestions that PT involves more than one water molecule. Here we point out that most of the studies that found a proton transfer through several water molecules have not involved a sufficient systematic search of the relevant reaction coordinates. This includes both energy minimization approaches as well as a recent metadynamics (MTD) simulation study. To illustrate the crucial need of exploring the potential surfaces reliably, rather than relying on automated approaches, we present here a very careful study of the free energy landscape along a 3D reaction coordinate (RC) exploring both the standard 2D RC, comprised of the attacking and leaving group reaction coordinates, as well as of the proton transfer (PT) coordinate. Our study points out that QM/MM minimization or MTD studies that concluded that the hydrolysis of phosphate monoesters involves a PT through several water molecules, have not explored carefully the single water (1W) path (that involves a direct PT form the attacking water molecule to the phosphate oxygen). Furthermore, we identified the most likely reason for the difficulty in finding the 1W path by QM/MM minimization methods, as well as by the current MTD simulations. We also discuss the problems with current studies that challenge the phosphate as a base mechanism and emphasize that all recent studies found associative/concerted paths (although many have not realized the meaning of their results). Finally, although we clearly do not have the last word about the 1W versus 2W paths we believe that we illustrated that the crucial mechanistic problems with alternative pathways should not be resolved by just running black box search approaches.

Exploring, refining, and validating the paradynamics QM/MM sampling.

The performance of the paradynamics (PD) reference potential approach in QM/MM calculations is examined. It is also clarified that, in contrast to some possible misunderstandings, this approach provides a rigorous strategy for QM/MM free energy calculations. In particular, the PD approach provides a gradual and controlled way of improving the evaluation of the free energy perturbation associated with moving from the EVB reference potential to the target QM/MM surface. This is achieved by moving from the linear response approximation to the full free energy perturbation approach in evaluating the free energy changes. We also present a systematic way of improving the reference potential by using Gaussian-based correction potentials along a reaction coordinate. In parallel, we review other recent adaptations of the reference potential approach, emphasizing and demonstrating the advantage of using the EVB potential as a reference potential, relative to semiempirical QM/MM molecular orbital potentials. We also compare the PD results to those obtained by direct calculations of the potentials of the mean force (PMF). Additionally, we propose a way of accelerating the PMF calculations by using Gaussian-based negative potentials along the reaction coordinate (which are also used in the PD refinement). Finally, we discuss performance of the PD and the metadynamics approaches in ab initio QM/MM calculations and emphasize the advantage of using the PD approach.

Paradynamics: an effective and reliable model for ab initio QM/MM free-energy calculations and related tasks.

Recent years have seen tremendous effort in the development of approaches with which to obtain quantum mechanics/molecular mechanics (QM/MM) free energies for reactions in the condensed phase. Nevertheless, there remain significant challenges to address, particularly, the high computational cost involved in performing proper configurational sampling and, in particular, in obtaining ab initio QM/MM (QM(ai)/MM) free-energy surfaces. One increasingly popular approach that seems to offer an ideal way to progress in this direction is the elegant metadynamics (MTD) approach. However, in the current work, we point out the subtle efficiency problems associated with this approach and illustrate that we have at hand what is arguably a more powerful approach. More specifically, we demonstrate the effectiveness of an updated version of our original idea of using a classical reference potential for QM(ai)/MM calculations [J. Phys. Chem. 1995, 99, 17516)], which we refer to as paradynamics (PD). This approach is based on the use of an empirical valence bond (EVB) reference potential, which is already similar to the real ab initio potential. The reference potential is fitted to the ab initio potential by an iterative and, to a great degree, automated refinement procedure. The corresponding free-energy profile is then constructed using the refined EVB potential, and the linear response approximation (LRA) is used to evaluate the QM(ai)/MM activation free-energy barrier. The automated refinement of the EVB surface (and thus the reduction of the difference between the reference and ab initio potentials) is a key factor in accelerating the convergence of the LRA approach. We apply our PD approach to a test reaction, namely, the S(N)2 reaction between a chloride ion and methyl chloride, and demonstrate that, at present, this approach is far more powerful and cost-effective than the metadynamics approach (at least in its current implementation). We also discuss the general features of the PD approach in terms of its ability to explore complex systems and clarify that it is not a specialized approach limited to only accelerating QM(ai)/MM calculations with proper sampling, but rather can be used in a wide variety of applications. In fact, we point out that the use of a reference (CG) potential coupled with its PD refinement, as well as our renormalization approach, provides very general and powerful strategies that can be used very effectively to explore any property that has been studied by the MTD approach.

Molecular thermodynamic modeling of the morphology transitions in a solution of a diblock copolymer containing a weak polyelectrolyte chain.

For a solution of the diblock copolymer composed of a hydrophobic block and a weak polyelectrolyte block, we obtain regions of stable aggregate morphologies in pH-solution salinity plane with the aid of the self-consistent field theory in the strong-segregation approximation. Lamellar, cylindrical, branched cylindrical, and spherical aggregates have been considered in the large interval of pH and salinity. The morphology stability maps are obtained to help control self-assembly of aggregates by variation of pH and salinity of the medium. In qualitative agreement with experiment, our calculations predict the coexistence of long wormlike micelles with branched and spherical micelles in transition zones. We compare the results of our calculations with available computer simulation and experimental data on micelles and brushes (planar and curved) formed by a diblock copolymer with one polyelectrolyte block. We show that for both weak and strong polyelectrolytes the agreement between the theory and experiment is satisfactory in most systems.